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Foreword 


G.  I.  Taylor  found  a solution  of  the  pas  dynamic 
eouaticns  which  describes  the  pressure  waves  produced  by  an 
explosion.  This  description  can  be  used  at  distances  from 
the  explosive  large  compared  to  the  explosive  dimensions, 
but  small  enough  that  the  shock  pressure  is  large  compared 
to  atmospheric  pressure.  In  fact  the  explosive  is  assumed 
to  be  a point,  and  the  atmospheric  pressure  is  neglected 
compared  to  the  shock  pressure,  in  obtaining  the  solution. 

For  these  reasons,  the  solution  is  called  either  the  "point 
blast"  or  "strong  shock"  solution.  The  same  solution  was 
also  obtained  by  J.  von  Neumann,  and  a solution  of  the  same 
tyre  was  found  for  explosions  in  water  by  H.  Primakoff. 
Similar  solutions  were  studied  by  G.  Guderley  in  Germany. 

The  two  opposing  restrictions  or.  the  range  of  validity 
of  the  Taylor  "point  blast"  solution  are  such  that,  for 
ordinary  explosives,  there  is  practically  no  range  in  which 
they  are  both  satisfied.  For  nuclear  explosives,  however, 
the  size  of  the  explosive  is  so  small  tnat  there  is  a range 
in  which  this  solution  is  useful.  This  has  been  demonstrated 
by  comparing  this  solution  with  experimental  results. 

In  order  to  extend  the  solution  to  greater  ranges, 
v;ho re  the  shock  is  weaker,  it  Is  necessary  to  tuKe  account 
of  the.  atmospheric  pressure  ahead  of  the  shock.  The  present 
report  by  Dr.  Morawetz  attempts  to  do  this  by  determining  a 
correction  to  the  "point  blast"  solution,  which  correction 
is  due  to  the  previously  neglected  atmospheric  pressure. 

A system  of  linear  equations  is  obtained  for  the  correction 
and  their  solution  is  analyzed.  Finally  for  the  case  of 
Primakoff’ s "point  blast"  in  water,  the  solution  is  given 
explicitly  and  graphs  of  various  quantities  are  given  for 
both,  the  original  o.nd.  the  improved  solutions. 

It  cannot  be  expected  that  the  range  of  validity  of 
the  "point  blast"  solution  can  be  extended  very  much  by  such 
a perturbation  method.  However  the  qualitative  and  quantita- 
tive nature  of  the  corrections  can  give  an  indication  of  the 
range  at  which  the  solution  begins  to  fail,  and  can  also 
indicate  the  manner  in  which  the  solution  changes  due  to 
atmospheric  pressure. 

Joseph  3.  Kc.llcr 
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1. 


PER  TITRUAT  TOUR  ABOUT  STRONG  SPHERICAL 


SHOCK  VaVES 


Introduction. 

It  has  been  shown  by  G.  I.  Taylor  and  others  that  there 
are  solutions  of  the  equations  of  spherical  flow  of  the  form 


u = X'1  £ Uc(r"Xt) 
c = X 1 T Co(r  Xt) 

D = x'2  ;4  ?0<p“Xt> 

where  \ is  any  constant  and  U . C , P satisfy  nonlinear 
ordinary  differential  equations.  These  solutions,  however, 
can  represent  a flow  behind  a shock  only  if  we  can  neglect 
the  pressure  ahead  of  the  shock.  Here  we  shall  study  the 
first  order  effects  of  this  pressure  and  the  modifications 
produced  on  the  original  flow. 

The  problem  of  finding  the  flow  behind  an  expanding 
shock  wave  of  constant  energy  is  reduced  to  solving  some 
ordinary  differential  equations  which  depend  only  on  y,  tne 
ratio  of  the  specific  heats.  In  the  case  of  y - 7j  these 
equations  can  be  solved  and  the  perturbations  in  the  flow 
quantities  expressed  explicitly. 

The  difference  between  the  flow  behind  a strong  shock 
and  a strong  detonation  can  also  bo  studied  in  the  same  way. 


lie  Equations  of  Notion . 

The  equations  for  spherically  symmetric  flow  are: 


u,  + uu  + — p =0 
t r p r 


(1) 


2u> 


Pt  + uPr  + P(ur  + -p) 


(rp"T)t  + u{pp“Y)r  = 0 , 

G.  I.  Taylor:  The  formation  of  a blast  wave  by  a very 

intense  explosion.  Proceedings  of  the  Royal  Society.  Series  A, 
Vol . 201,  March  1950,  n.  159. 
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where  u is  the  radial  velocity,  p is  the  pressure,  p is  the 
density  and  y is  ratio  of  the  specific  heats  of  the  gas, 
while  r is  the  distance  from  the  origin  end  t is  time. 

If  we  introduce  the  variables 


(2) 

(3) 


>1  = r'xt 


U = >tr~"u 


J P 


C = Xtr  c — X 
2 2 -2 

P = Xtr  *p 


where  X is  any  positive  constant,  we  obtain,  using 

S3  = -x  *£.  + 4* 

In  ( 1 ) . 

yi  P 

HU^(l-U)+  UtUt  - y”1C2  -~3-  + U(aU-l)  +2aY_1C2  = 0 

tpt  fcct 

(5)  - +— jrMl-U)  - 2 — -Ml-U)  +U-p-  -2U  + 3aU  = 0 


^-(  1-0)  - ^u-u)  + o ^ - —f  0 


t + t§t(1“gTJ) 


= o 


where  a = A-^. 


Equations  (5)  may  be  rewritten  as 

-1  2 tPt 

= A - (l-U)tTJt  - y XC 


(6)  D 


c 


Y-l  J.TT  _ D tCt 

"T"  tut  T=TT  “ 


(y-l)  c2  tpt 

T^u  — 


V\P*  tpt 

D = E - ytUt  - (1-U)  -pi- 


whore 


r*~i  'WIO.- 


- ' c-  . . 
*••«.*  •?.*  •< 
. . . . 


TV's  !*&■' 
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A = ( 1-U) (l-aU)U  ~ ( 3aU-2Y"1(l-a) )C2 

B = ( 1-U)  (1-aTT)  -(Y-l)[(l-aU)  - ^(l-o)JU-  [G+Y'1(l-a)(l-TT)"1]c2 
E = ( 1-U) i 2( 1+y) ( 1-aU)  - 3y]  + YU-aU)  - 2aC^ 

D=(l-U)2-C^ 

There  are  two  identities  relating  A,  B,  D and  E, 


(8) 


-A  + B(l-TJ)  + D(3aU-^§Y  (1-aU))  = 0 
-A  + y_1  E(l-U)  + D(3aU-2  (1-aU))  = 0 


The  special  solutions  for  spherical  flow  known  as 
spherical  "progressing"  waves  are  solutions  of  (6)  that  are 
independent  of  t.  They  represent  states  behind  or  ahead  of 
shocks  and  detonations  if  the  shock  or  detonation  front  is 
given  by  Y\  = constant.  We  represent  any  such  flow  by  U0(t^)  , 
CQ(  i\)  t p0(^)  • Then  from  (6),  we  find  the  following 
nonlinear  ordinary  differential  equations  for  UQ,  CD  and  Pq, 


(9) 


n 


0n 


f?- 


E. 


where  Aq  = A(U0,C0)  etc. 

Prom  the  first  two  equations  of  (6)  we  obtain  an 
equation  for  CQ  as  a function  of  UQ, 


i 


P 


o 


dC0  . 


3Cn 
0 o 


The  solutions  of  this  equation  then  yield, 
and  as  functions  of  UQ. 


through  ( 6) , 


We  shall  consider  here  flows  that  arc  approximately 
progressing  wave3,  that  is  flows  given  by 


(10)  U = Uo(y^)  + eU1(^,t)  + 0(e2) 

C - C0(yl)(l+  eC1(  y^>t)  + 0( e2 ) ) 

P = Pc{  n.)d+  eP1(  n,t)  + 0(e2)) 


where  e is  some  small  parameter. 

Substituting  (10)  in  (6)  and  equating  first  order  terms 
in  e we  find 


ui+  (b)c  coci  -D^U-VWit- 

ui  + (l)°  coci  - *ir  D;lwit 

. .Vifit.  eL„-iw 

Ul+  (It  CoCl-  U-plD^tPit  . 


D 


u=uolyi),  c=c0(»v) 


Conditions  at  a Shock . 

The  boundary  conditions  on  the  flow  quantities  at  a 
shock  are  given  by 


p(u-z)  = - p^z 
p(u-Z;  + p = 


(12) 


2 

p-i  z + p 


5. 


where  z Is  the  velocity  of  the  shock  front,  and  p^,  p-^  are 
the  density  and  pressure  ahead  of  the  front. 

In  terms  of  the  variables  U,  c,  P of  (3)  equations  (12) 


become 


o(U-Z)  — - p^Z 


/ t ...  „,2iD  _2  .,2  „2a.  2-2a 

(13)  pin-z)  +p  = p1z  + x t Pj 


2 , T _ * 2 , - 2._2  2 2 . ..  2.  2, 2 2a,  2-2a 

p,  (U-Z)  + ( 1-p  )C  = p,  Z + (1-p.  )c^a  t 


where 


Z :=  Xtr  z 


(15)  p2  = (y-iHy+D"1 


2 _ -1 

C1  ~ YPiPi 


From  (13)  we  see  that  a spherical  progressing  wave  can 
represent  the  state  behind  a shock  for  X f 1 only  if  c-^  = p1=0, 
that  Is,  only  if  the  shock  is  infinitely  strong.  In  this 
case  the  position  of  the  front  is  given  by 

un  i=h0 

where  Hq  Is  some  constant  and  thus  from  (14) 

(18)  Z = Xtr"1  = 1 

We  assume  now  that  the  deviation  in  the  shock  path  is 
small,  and  set 

(19)  ft  = HQ(l  + eH1(t))  on  the  3hock 
Then  from  ( 14)  we  find,  to  first  order  in  e 

(20)  Z = 1 - etH,*. 


Substituting  (2o)  in  (13)  and  setting 


(21) 


6 


wc  obtain 


Vi  = (l-p.2)(l  - etKlt  -er2X"2) 


( £2 } 


C2  - u2(  1+1X^ ' ( 3-  - 2eoi-H 


lr 


l-2u^  *2^2a, 

+ e -M  giy-  t.  ) 


^(l+l/) 


p = pT(l-ti2)(l-  2etHit  -e  -±~z  t2~da) 


1+|A 


for  H = Hc(l+  eH1(t) ) . 

On  the  other  hand,  expanding  TT(H0(l+eH1(  t)  ,t,e)  etc.  in 


powers  of  e we  find 


(£3)  t:(h  (l+e:i.,(t)),t,e) 


= uo(Hn(1+eHl  ( t) )”)  + eu1(H0(l+eHi(  t) ) ,t)  + 0(  e ) 


= VHo}  +enl(I'o»t)  + eHlHo  _ Hl(t)+0(e2) 

n-H0 


'A 


= V^J  + eU^ti  + e (ir 


H^( t ) + C( e ) 


C(H0(l+eH1(t)),t,e) 

= C (H0)  jl  + eCn(II0,t)  + e J ^(t)  + C(e2) 

L ‘ v 0/n=Ho 

P(H0(l+eH1(t) ),t,e) 


= r0(H0l 


i+  Ep,;s0,t)  +£(d;) 

' * \A 

Comparing  (23)  with  (22)  wc  have 


H1( t)  + 0( e2) 


n=H, 


TT0(H0)  = 1 - 


(2k) 


P p ? . 

C$(K_)  = n“(l  + M-  ) 


Ji 


di 


»M^ 


It  is  then  clear  from  (9)  that  UQ,  CQ  and  FQ/p-j_  are  all 
functions  of  *VH0  which  depend  only  on  y* 

From  the  first  order  terms  in  e,  we  find 


Ul(Ho,t)  = 


(25)  C1(H0,t)= 


H^( t)  - (l-|i2)tHlt(t)  - ( 1— (i2) t2-2a 

V °\=H0 

- (jr ) hi‘  ‘>-  tHit'  + ■ 1:^2.  t2‘2a 
Vo/  2(4,  ( 1+^1  ) 


pllH0*t)  ~ ' (l5^  ) Hl(t)  “ 2tHlt(t)  “ ^ t2_2a 


The  perturbation  H-^(t)  in  the  position  of  the  shock  can 
be  eliminated  and  (25)  reduced  to  two  conditions  on  the  flow 
quantities,  namely  for  Y[-  HQ, 

p p 

ttl^l  + (^1^1  Yi^t q ~ 

(26) 

a2Ul  + ^2C1  + y2P1  + a21tUlt  + ^21tClt  = 62fc 


where  a-^,  etc.  are  all  constants. 

To  find  the  flow  behind  a shock  we  now  have  to  solve 
the  system  of  hyperbolic  differential  aquations  (11)  of  third 
order  and  two  conditions  (26)  on  a space-like  line,  the 
shock  path.  One  more  condition  must  be  prescribed.  For 
example,  we  might  prescribe  the  velocity  cf  a particle 
corresponding  to  a given  piston  motion  which  does  not  differ 
much  from  that  which  maintains  a progressing  wave.  This 
would  involve  extensive  computations  although  it  can  be 
reduced  to  a second  order  problem. 

Alternately  and  more  naturally,  we  may  prescribe  the 
total  energy  contained  in  the  shook  wave,  that  is,  the 
energy  imparted  at  t = C.  This  Is  essentially  equivalent 
to  prescribing  one  condition  on  the  t-axis,  namely  that  no 


t;  •- J • 

WrUigWuii  laa  H i 


■ A -*•-  n 

y.  ■ t - • * • 
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energy  is  addei.  Prom  the  general  theory  of  hyperbolic 

equations  we  will  have  trie  correct  number  of  conditions  to 

determine  the  problem.  In  fact  this  problem  can  be  solved 

in  terms  of  the  solutions  of  ordinary  differential  equations. 
In  the  case  of  a shock  wavo  in  water-  it  can  bo  solved 
explicitly . 


Shock  Wave  of  Constant.  Energy . 

First  we  note  that  equations  (11)  with  tho  boundary 
conditions  (2w)  havo  the  special  solutions, 


TJi  = t2_2a(  )/n{  y{)  + h y.12<  n>> 
C;L  = t2"2a(X.2i(  + h X2?(  n.)) 

?!  = t2_2a(x3i(n>  + hX32<  n.)) 


(27) 


H,  = ht 

X 


2-2a 


where  h is  an  arbitrary  constant  and  Xn»  X21*  X3i  *3 
solution  of 

4“  ■ (b L Xi  + (c)°  coX2 


<“)  - 01  * + (§)"  co^  - 


v " /U 


2 1-a 


ci;1 


l-TI  C 
o o 


lT5f-M^  + MCoX2-2r(x-.)l»S1Xi 


.-1, 


" 2(l-a)D0A(l-U0)  V3 


which  satisfies  the  initial  conditions 


- f3pWt.  . 


-•  V*  '•  Vv-  w-  ‘ 


5536#--  „ 


(29) 


X21<Ho> 


- 2, 
-l-L-li  ) 


L-2u^ 


2^‘"(  l+p.^7 
2 

_IE 


l+p' 


9. 


and  Xi2»  ^-22*  X32  ~s  solution  of  (28)  which  satisfies 

the  initial  conditions 


(30) 


Xl2<Ho>  = * (c;)  - 2<l-|/)(l-a) 

X22^Io^  " “ 


X32(ho>  = - (ir)  -W1-*)  • 

° VHo 


A method  for  reducing  the  third  order  system  (28)  to 
one  of  second  order  is  contained  in  the  Appendix. 

It  turns  out  that  the  solution  (27)  is  Just  the  solution 
which  satisfies  the  condition  of  constant  energy  if  h is 
chosen  appropriately. 

Let  the  total  energy  contained  in  a spherical  shod: 
wave  at  any  time  be  E(t)  where 

2(t7i  2 p-pi\  2 

(31)  E(t)  = kn  j pu  + -y=T)  r <*r 


Hero  R(t)  is  the  position  of  the  shock.  In  terms  of  the 
variables  (2)  and  (3)>  (3D  becomes 

(32)  E(t)  » -I4.no3  J*  l2j,  + P tSa-2  av^ 

where  1^  = H(t)  represents  the  shock. 


In  the  ease  of  an  infinitely  strong  shock,  where  = 0, 
E(t)  will  be  constant  only  if  n = 2/5.  Then 

(33)  E=E  0 = ^!f  f ( z ~pr  + v=t)  Po  1~3  d *1  • 

% V Co  1 

If  we  now  consider  that  second  order  terms  in  can  be 
neglected  and  use  the  approximations  (10)  and  (19)  we  obtain 

<3«  E = - if  J°  (2  + 7^1)  n"3  dn 

00  o 

' 2 1 

f -i  YU„  -1  \ _ •> 

" ta£  eHoHi  f £ + po  n 

L 0 irl=H0 

” "2fif  f - (f  ^1)  popi+  ~r(urTJoci)i>r3  dri+  °^e2 

't  o o J 


= + nrs r-r  n ' ' z ' ' . 

0 3 ( Y“ 1 ■ 0 

From  the  first  order  terms  in  e using,  from  (21), 


1 „-6/5  ,6/5 


« - X2  4 H 2ca  - # Wip/  hVS 


we  obtain 


, f /,  rD?  \ 1 ,%  f / u?  , \ 

551 


!tl=H0  -co 

yu^p 


o o,. 


Ul‘Uocl)f  V 


pl  .6/5  tt-2 

TiFTT  - Ho 


or,  using  ( 2lp)  , 


11 


(36)  H,  SU=£  ♦ £ { (?  | * £)  P, 

Yur,.  . fix  y*3  /w  ^ _ 5 1 *.6/5 

+ °,VCo0i)S--^)  <>(£)--  frTFTT1 

A 

W / 

Wo  can  new  show  that  (27)  is  just  tho  solution  of  (xx) 
and  (25)  which  satisfies  (36)  if  h is  chosen  appropriately. 

Substituting  (27)  in  (36)  using  r ^t  = H0  to  first  order 
on  tho  shock  yields 

1 r 


ht 


6/5 


(*!•*)** 
,-21/5 


~f(  X12-uo  X22>J  ^ (^)’21/5d  (O] + *6/#  ^ + ^r) 


yU 

+ — £( 

Co 

or 


X iruo  X21>)  57  (^)‘21/S  d (0  = - i rT^n 


.6/5 


(37) 
where 
A 


A-jh  + A2  = A^ 


woo  L v wo 


yU 

c 


O . w 

' w 

v A 12 

O 


..  ■./  A ?o  m t21/5  d ('il'i 
•uo'*2'/  V*J 


n f/  U2 


yU, 


(33)  a2  = J X31+-t?(Xxi-d°  X2i> 

a3  = - I TTrTT 


12. 


Equation  (3?)  can  always  be  solved  for  h provided  A-j  ^ 0. 

Thus  the  first  order  increment  in  the  energy  vanishes  and 
we  have  a complete  solution  of  the  problem  in  terms  of  Xu 
etc . 

From  (7)»  (9)  and  (24)  we  see  that  and  CQ  and  PQ/ p-^ 
are  functions  of  Y(/  H0  only  and  independent  of  while 
P/p^  is  a function  of  • 

Thus  the  coefficients  of  the  equations  (28)  are  functions 
of  *l/Hn  and  hence  by  (28),  (29)  and  (30)  the  functions  Xu* 

*21’  %3i»  Xi2»  XzV  X&  dePend  on  *tA* o only. 

Thus  and  Ag  are  constants  depending  only  on  the 
solutions  of  fixed  differential  equations  with  fixed  initial 
conditions,  depending  only  on  y. 

Substituting  (27)  in  (10)  and  then  in  (3)  we  find 


It  i'Io+Et6/s‘Xii+hXu'> 

| t-3/S  n-2/S„o  + 1 et3/5  n-2/5<  X u + h X12) 

f t-3hn-2/5Co+|  Et3/S  21* 


I.  4 /a  _ !i  /c*  I. 


»f"  "po  + % E *1  H/  "f  o(  x 31  * h X32) 


Note  that  for  fixed  the  perturbations  in  u and  c 
increase  with  r while  the  perturbation  in  p is  constant. 

Now  the  pressure  dies  out  behind  the  shock  and  the 
maximum  pressure  occurs  at  the  shockj  thus  from  (22)  and  (3) 
we  find 


(J£  h+-^)(l-n2)P1 

^ 1+u  1 


In  other  word 3,  the  maximum  pressure  is  the  pressure  of  the 
Taylor  point  blast  wave  plus  a constant. 


(4-0) 


13. 


Shockwave  of  Constant  Energy  In  Water . 

In  the  case  of  an  explosion  in  water  the  functions 

etc.  can  be  calculated  explicitly-. 

Here  we  have  F 3 -1)  where  A = 3000  atmospheres 

Po 

and  y = 7»  Then  we  replace  the  last  equation  of  (3)  by 
P = X2t2r~2(  p + A) . In  this  case, 

(39)  UQ  = 1-n2  , C2  = iAi+m2)  , P0=  Pl(l-ti2)  V2A 

is  the  solution  of  (9)  and  (24)  . Then  An  = 0,  B = 0 while 
Eq  and  Dq  are  constants. 

The  differential  equations  (28)  then  have  constant 
coefficients  in  the  lowest  order  terms  and  the  solutions  are 
powers  of  . Satisfying  equations  (29)  and  (30)  \jg  find. 
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This  special  solution  is  due  to  H.  Primakoff. 
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In  Figures  1,  2 and  3 wo  have  plotted  the  path  of  the 
shock  as  a function  of  time  and  the  maximum  pressure  as  a 
function  of  time  and  distance. 


State  Behind  a Detonation . 
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include  tho  chemical  energy  of  the  detonation.  Wo  then  have 
for  the  conditions  across  the  front 


(1+2) 


p(u-z)  = - p.^ 


2 2 
p(u-z)  + p = p-^z  + p-j^ 


T-l^V"r  + £(U-z)&  + E — ~rr~  + " ± , \ 

(r-1)  2V  2 p^i,  y-^-1) 


>2  . “ 


YP-i 


+ Kn 


whore  E and  are  the  energy  of  formation  per  unit  mass  of 
the  burnt  and  unburnt  material  respectively,  and  y^  is  the 
ratio  of  the  specific  heats  in  the  unburnt  gas.  In  torms  of 
the  variables  (3)  and  (lip)  the  last  condition  of  (1+2)  becomes 
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If  we  perturb  about  a strong  detonation,  i.e.  setting 

,2  2 ^ — 

it-r(  1-li ) + E - E,  = 0 we  can  find  the  undisturbed  flow  as  a 

c ' ^o  o 1 

spherical  wave  and  the  perturbed  flow  satisfies  (25)  with 

2\«»2 

different  constants  for  the  coefficients  of  r . 

In  this  case  there  are  again  special  solutions  of  the 
form  (27)  where  now  the  initial  conditions  (29)  must  be 
adjusted  appropriately.  However,  in  this  case  we  need  one 
more  condition  on  the  flow  and  the  special  solution  will  not 
in  general  satisfy  it. 

For  y = 2 we  have  U (H„)  = CQ(H0)  or  u = c at  the  shock. 
The  unperturbed  detonation  ia  then  a Chapman- Jouguet 
detonation.  If  the  perturbed  flow  is  also  behind  a Chapman- 
Jouguet  detonation  we  obtain  a relation  between  TJ1(H0),  C-,(HQ) 
and  H^(r).  This  relation  can  bo  satisfied  by  solutions  of 
the  form  (27)  for  an  appropriate  choice  of  h. 


Appendix. 

In  this  Appendix  we  will  show  how  equations  (28)  can  be 
reduced  to  one  second  order  equation.  In  general,  equations 
(6)  and  the  shock  conditions  (25)  can  be  reduced  to  a second 
order  equation  involving  the  unknown  shock  function  H^(x). 
From  equations  (8)  wo  see  that 

r^=  0 
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and  from  (11)  then, 
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(A.2) 
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If  wo  introduce  new  depondont  variables, 

(A. 3)  «x(  VI)  - - Xi(  AJ/V1’^) 

V\)  = y=T  X2(  *l>  + *1(  *l)Uo 

«3(n)  = y"1  X3<*i>  + ^i(n.)u0 

equations  (A.2)  reduce  to, 
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u~  " {'3  + y=T 


U-l)  + 2(l-a))£,  - (2-2a)§; 


1-V„  d*. 


u. 


= (-3  + 2Y‘1U-1>  + 2(l-a))£-  - TJ~1(2-2a)% 


By  introducing  now  dependent  variables  we  can  reduce  (A.ij.) 
to  a much  simpler  form,  Wc  set 

(a. 5)  y = f — 

i i_vo  -v 


where  H#  is  any  fixed  value  of  Thon  (A.ij.)  roduces  to 

ae,  C2 

(A. 6)  D, 
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'o  K0^1  + V2  + 2(i“a)  dr  ^3  55  0 
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where  KQ  and  LQ  are  functions  of  y, 
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and  jig  and  are  constants. 
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(A. 8}  sin  - 3 “ 2v  1?v  + 2v  *"  - 2( 1-a) 

^3  = 3 - X + —j  - 2(l-a) 

•The  solutions  of  the  third  order  system  (A. 6)  con  also 
be  expressed  in  terms  of  the  solutions  of  a second  order 
differential  eqviation.  Wo  introduce  tho  function  G where 


(A. 9) 


= € 


1 


(A. 10) 


G(YC)  = 0 


whero  Y is  some  fixed  value  of  y. 
o 

Then  from  the  last  equations  (A. 6)  wc  obtain 


( A . 11 ) 


£2(y)  ~ ^2(Yo*  + ^2G^y- 


i 

! 


(A. 12)  €3(y)  = C3(Yo)  + ^G(y) 

i 

Substituting  (A. 11)  and  (A. 12)  in  tho  first  equation  of  (A. 6) 
yields 

2 

(A. 13)  D0  + Ko  ||  + g(h2L0  + (2-2a)^3  } 

ay  ^ u~  / 

2 i 

= -W-2i1-i7«3'V  • ■ 

o 

Now  an(^  ?3  are  prescribed  for  y = iQ  by  (29)  or  (30). 

Thus  we  have  a differential  equation  (A*13)  for  G and  two  - 

boundary  conditions  (A»9)  and  (A.  10)  at  y = 

■--4 
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